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Abstract 

An integrable quantum spin ladder based on the SU(4) symmetry algebra with 
boundary defects is studied in the framework of boundary integrability. Five non- 
trivial solutions of the reflection equations lead to different boundary impurities. In 
each case the energy spectrum is determined using the quantum inverse scattering 
method. The thermodynamic properties are investigated by means of the thermo- 
dynamic Bethe ansatz. In particular, the susceptibility and the magnetization of 
the model in the vicinity of the critical points are derived along with differing mag- 
netic properites for antiferromagnetic and ferromagnetic impurity couplings at the 
edges. The results are applicable to the strong coupling ladder compounds, such as 
Cu 2 (C 5 H 12 N 2 ) 2 Cl 4 . 

Key words: Yang-Baxter equation; reflection equations; thermodynamic Bethe 
ansatz; impurity effects, quantum spin ladders 
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1 Introduction 

Research on quantum spin ladders continues to attract considerable attention 
from both theoretical and experimental points of view due to their relevance 
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to a large number of low-dimensional materials, such as particular cuprates 
and organic compounds [1,2,3,4,5,6], among others. Initially, most of the the- 
oretical results concerning ladder systems were obtained from the standard 
Heisenberg ladder. Subsequently, other models with generalized interactions 
have been proposed. In this context, Nersesyan and Tsvelik [7] introduced 
a spin ladder model incorporating a biquadratic spin exchange interaction 
term, which, when sufficiently strong, exhibits new dimerized phases [8]. Var- 
ious ladder models have been developed by an extension of the symmetry 
algebra [9,10,11,12,13]. A special case of the Nersesyan-Tsvelik model [7] was 
proposed later by Wang [12]. This model, based on the S77(4) symmetry al- 
gebra, is exactly solvable by Bethe ansatz methods and exhibits a spin gap 
in the spectrum of elementary triplet excitations, a necessary condition for 
superconductivity to occur under hole doping. In addition, it was recently ob- 
served [14] that this model can be used to describe some physical properties 
of different types of two-leg ladder compounds, such as Cu 2 (C5H 12 N 2 ) 2 Cl 4 [3], 
(C 5 H 12 N) 2 CuBr 4 [4], (5IAP) 2 CuBr 4 -2H 2 [5] and KCuCl 3 , TlCuCl 3 [6]. In the 
absence of a magnetic field the model exhibits three quantum phases, while in 
the presence of a strong magnetic field there is a gapped phase in the regime 
H < H cl , a fully polarized gapped phase for H > H c2 and a Luttinger liq- 
uid magnetic phase in the regime H cl < H < H c2 . This observation suggests 
that the physical properties of the ladder compounds can be accessed via the 
well-established knowledge of integrable systems. 

On the other hand, the effect of boundary impurities and defects also plays an 
important role in quasi-one dimensional systems. An integrable S77(4) spin 
ladder model with a boundary defect has been proposed and investigated 
recently through Bethe ansatz methods [15]. This model, however, is just a 
particular case of a more general family of exactly solvable ladder models based 
on the 577(4) symmetry algebra that can be constructed from more general 
types of bounday conditions. Basically, by this strategy, a set of equations 
to deal with the boundaries, called reflection equations (RE) are introduced 
[16,17]. The solutions of these equations [18,19,20], referred to as TT-matrices, 
in turn introduce boundary interactions into the Hamiltonian of the system, 
in such a manner that integrability is preserved. The boundary interaction 
terms in spin ladder models may be realized by impurity doping at the ends 
of the ladder. Impurity doping in a spin ladder system with a spin gap has 
been performed [21]. Substantial change in macroscopic properties such as 
enhancement in spin correlations and magnetic susceptibilities are observed 
in the low impurity concentration region. The boundary impurity doping may 
change the critical behaviour at the boundaries of the ladder systems. 

The purpose of this paper is to present a complete family of integrable spin lad- 
der systems based on the SU (4) symmetry algebra with boundary impurities 
in a systematic way. An analytic analysis of the thermodynamic properties of 
these models is then performed by means of the thermodynamic Bethe ansatz 
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(TBA) method and, in particular, the effect of these impurities on the free 
energy, the susceptibility and the magnetization is discussed. So far, the re- 
sults obtained provide a clear interpretation of the impurity effects in the low 
temperature regime of an integrable open spin ladder system. 

The paper is organized as follows. In section 2, we present the SU(4) solu- 
tion of the YBE and solve the corresponding RE. Furthermore, we give the 
explicit expressions for the Hamiltonian with different types of boundary de- 
fects. Section 3 is devoted to the derivation of the Bethe-ansatz solution by 
means of the Quantum Inverse Scattering Method (QISM). The reader more 
interested in the physics of the model may choose to skip section 3. In section 
4, the ground state properties, the quantum phase diagram and the boundary 
impurity effects are studied via the TBA. A summary and discussion of our 
main results is given in section 5. 



2 The integrable spin ladder model with boundary impurities 

Let us begin by introducing the integrable spin ladder model based on the 
SU (4) symmetry with boundary fields, 

H = ^H leg + J ± f: ■ % + H™ + H%>, (1) 
7 j=1 

where the leg part consists of Heisenberg exchange and four-spin interaction 
terms 

^ieg = + ' ^J'+i + T i ' Tj+i + • Sj+i ■ Tj ■ T j+1 j . (2) 



The left (right) boundary terms H[ m ' (H L ') depend on arbitrary parameters 
U± and are given by 



H[ m) = { 



-U-^-fi - \U-, for m = 1 

-U-(\-Sl)(\-T?) + \U-, for m = 2 

f/_( 1 S 1 -f 1 -(|- 1 S 1 2 )(|-Tf)) + if/_,for m = 3 , 

[/_ (S, ■ f x - (I + S\){\ + T*)) + \U_, for m = 4 



-2USITZ + W-, 



for m = 5 



(3) 
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Fig. 1. The su(4) spin ladder with boundary impurities. U± are the boundary im- 
purity coupling constants. J» and J± are the intrachain and interchain couplings. 



H 



(0 



-U+S L ■ T L - \U+, for 1 = 1 

- U+ (\-Si){\-Ti) + \U + , for 1 = 2 

U + (S L -f L -{\-Sl){\-Tl)) + \U + ,ior 1 = 3 

U + (S L -f L -{\ + Si){\ + Ti))+\U + Jov l = A 

-2U + S'iTl + \U+, for 1 = 5 



(4) 



— * — » ^ 
In the above Sj and Tj are the standard spin— ± operators acting on site j 

of the upper and lower legs, respectively, Jn and Jj_ are the intrachain and 

interchain couplings (see figure 1) and L is the length of the ladder. It is 

worth noticing that 7 is a rescaling constant which can be used to minimize 

the biquadratic term such that the quantum phase of the model (1) lies in the 

same Haldane spin liquid phase as that of the conventional spin ladder (see 

section 4). 

Notice that the boundary terms corresponding to the first case (m = 1, I = 1) 
in (3) and (4) act as Heisenberg-type rung couplings, whereas in the second 
case (m = 2, I = 2) they act as a ^-component spin interaction with boundary 
magnetic fields. In the third (m = 3, / = 3) and fourth (m = 4, / = 4) cases, 
they act as a combination of Heisenberg-type rung coupling and z-component 
spin interaction with boundary magnetic fields. In the last case (m = 5, I = 5) 
only the z-component spin interaction terms survive. The Hamiltonian (1) 
thus contains five different types of boundary rung interactions at each edge 
of the ladder realizing different impurity dopings. This leads to twenty five 
possible choices of boundary impurities. The rung interaction in the bulk was, 
as usual, introduced by the chemical potential terms given by —J±Y1 
in the canonical basis 
denoted by 



l n 
j=l e j 

Cj. The rung states split into a singlet and a triplet 



H> = ^(IU> 



IT)), |2) = |TT), 
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|3> = ^(|U> + UT», |4> = |U>, 



(5) 



respectively. The leg interaction part of the Hamiltonian (2) does not change 
under the basis transformation (5). However, the bulk rung interaction part 
and the boundary rung interaction terms alter with respect to the choice of 
the order of singlet and triplet in the basis (5). 

In order to derive this model let us begin by recalling the SU (4) i?-matrix 



R 12 (u) 
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(6) 



where the elements are given by 

Wi — u + 1, w 2 = u, w 3 = 1. 



(7) 



The quantum i?-matrix (6) satisfies the Yang-Baxter equation (YBE) 
R 12 {u - v)R 13 (u)R 23 {v) = R 23 (v)R 13 (u)R 12 (u - v), 



(8) 



guaranteeing the integrability of the model with periodic BC. This i?-matrix 
enjoys the properties 



Ri 2 (u)R 2 i(-u) = 1 - u\ R\ 2 t2 (u) = R 12 (u), 



(9) 



where superscript t denotes the transposition in the space with index a. 
For other types of boundary conditions, the YBE will still account for the 
integrability of the bulk part of the model, but the boundary terms have to 
be chosen appropriately in order to preserve the integrability. In particular, 
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the left and right reflection matrices, K_ and K + , respectively, are required 
to satisfy the RE's [17,19] 



Ri 2 (u - v) K- (u) R 2 i(u + v) K- (v) 
= K- (v) R 12 (u + v) K- (u) R 2 i(u - v), 
R 2 i 2 (v - u) K%{u) R 12 {-u - v) K%{v) 
= K% (v) Rl\ t2 (-u - v) K%(u) R% t2 (y-u). 



(10) 



(11) 



In the above we have introduced the object R which may be determined by 
the relations 



R%(-u)Rl\(u) = 1, Rl\(-u)R%(u) = 1, 



(12) 



and we have used the conventional notation 
X=X®l V2 , X=I Vl ®X, 



(13) 



where I v denotes the identity operator on V and, as usual, i? 2 i = P • Rn • P, 
with P being the permutation operator. After a lengthy calculation we find 
the possible solutions of the RE's for the diagonal if±-matrices (see also Ref. 
[18]) 



K±(u) 



(K\±(u) \ 

K2±(u) 

K3±(u) 

\ KA±{u)J 



(14) 



The solutions for K_, corresponding to the left boundary are: 



Case 1 

Kl_{u) = « + K2_{u) = K3_(u) = KA_(u) = -u + 



(15) 



Case 2 

Kl-(u) = K2^{u) = K3-(u) = u + KA_{u) = -u + 



(16) 



Case 3 

Kl_(u) = K4_(u) = — u + , K2_(u) = K3_(u) = u + 



(17) 
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Case 4 



Kl-(u) = K2_(u) = + K3_(u) = KA_{u) = u + 



Case 5 

= #3_(u) = u + #2_(u) = KA_(u) = -u + f_. (19) 

On the other hand, the solutions for if + , corresponding to the right boundary, 
are: 



Case 1 



Case 2 



if 1. 
KA 



+ 



Case 3 

ifl H 
if2 H 

Case 4 

ifl H 
if3 H 

Case 5 
K2 A 



«) = -« + £+ + 3, 

u) = K3 + (u) = KA + (u) =«-£+ + !; 



«) = A"2+(«) = if3+(w) 
u) = u - + 3; 



u) = KA + (u) 
u) = K3ju) 



u - + 2, 
-« - - 2; 



«) = J K'2 + (u) = u-£ + + 2, 
u) = KA + (u) = -«-£+ -2; 



«) = J K'3 + («) = -«-e + -2, 
w) = iT4 + (w) = + 2. 



(20) 



(21) 



(22) 



(23) 



(24) 



In the above £± = are free parameters related to the left — (right +) 

boundary coupling LL (£/+) respectively. The boundary pairs ifi" 1 ^), K^{u), 
I, m — 1, . . . , 5, lead to the boundary pair terms iji m ) and if^P in (3) and (4), 
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whose combination leads to twenty five possible choices of boundary impuri- 
ties. A special choice of the boundary term was already investigated in 
[15]. The symmetries enjoyed by the i?-matrix (6) and the RE's (10) and (11) 
constitute the necessary ingredients for the integrability of the model with 
boundary impurities, due to the fact that the double-row transfer matrix of 
the system 



t(u) = tr [^ + (m)T(m)K_(m)T- 1 (-m) 



(25) 



commutes for different values of the spectral parameter. Here T(u) denotes 
the monodromy matrix given by 

T{u) = R 0yL {u)Ro, L -i{u) ■ ■■ Ro, 2 {u)Ro,i (u) (26) 



and T 1 its inverse. The Hamiltonian (1) associated with the quantum R- 
matrix (6) is related to the double-row transfer matrix (25) by 

H = In r{u) | u=0 -J±J2 e f + const - (27) 



Here 

s^-'i- = + ^<°> A -<°> + 2 ^SP ' (28) 



where the prime denotes the derivative with respect to the spectral param- 
eter. The relation (27) clearly indicates the identification U± = between 
the boundary impurity couplings U± and the free parameters £± of the bound- 
ary scattering matrices. So far, we have completed the first step towards the 
solution of the model with boundary impurities. Next we proceed with the 
diagonalization of the transfer matrix (25) by means of the open algebraic 
Bethe ansatz [24,25]. 



3 The algebraic Bethe-ansatz approach 



3.1 First-level nesting structure 



In order to find the spectrum of our Hamiltonian with boundary defects, 
we first need to solve the eigenvalue problem of the transfer matrix, namely 
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t$ = A$. As usual, the transfer matrix (25) can be written in the form 



r(A) =tr K + (u)f_(u) 



(29) 



where T_ (w) is the double-row monodromy matrix defined by 

f_(u) = T{u)K_{u)T-\-u). (30) 



One can verify that T_(u) also satisfies the RE (10). Following the notation 
used in Refs. [22,23,24,25], we label the elements of the monodromy matrix 
T{u) by 



T(u) = 
and further 
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5 2 ( 
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D u 
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£>13 
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D 21 
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£>22 


(«) 


£>23 


:«) 






D 31 
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^32 


(«) 


£>33 


») 






«) 
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:«) 


S 2 ( 


«) 


5 3 ( 


u) \ 


)= 


Ci 


:«) 


Du 


(«) 


£>12 


(«) 


£>13 


:«) 






:«) 


D 21 


(«) 


£>22 


(«) 


£>23 


:«) 




V<?3 


:«) 


D 31 


(«) 


£>32 


(«) 


£>33 


:«)/ 



(31) 



(32) 



(33) 



According to the first level Bethe ansatz, the eigenvectors |$) of the transfer 
matrix can be written as 



\^)=B ll {u 1 )---R lN {u N )\ct ) )F^ 



'IN 



(34) 



where the summation is taken on the repeated indicies. The coefficients with 
indices % n = 1,2, 3, n = 1, • • • , N will be determined later by the second level 
Bethe ansatz. The first level pseudovacuum \<f>) is chosen as the standard 
ferromagnetic state 



|0>i 



•® |0>i <g> - - • <g) |0)] 



(35) 



where |0>* = (1,0,0,0)* acts as a highest- weight vector. This state corresponds 
to the product of the rung singlet state in the basis (5). Different choices 
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of the order of the basis (5) will change the eigenvalues of the Hamiltonian 
which facilitates the analysis of the ground state in different regions. From 
the structure of the i?-matrix (6), one can deduce that the operators B^u) 
and Bi(u) (i = 1, 2, 3) act on the reference state as creation operators creating 
particles with pseudo-momenta u and —u, respectively. The operators Ci{u) 
(i = 1, • • • , 3) behave as annihilation fields. Furthermore, using an invariant 
version of the Yang-Baxter algebra, 

T-\-u)R 12 (2u) T(u) =T(u)R 12 (2u) T'^-u), (36) 



we obtain, apart from an overall factor Q(u) = Kl_(u)Kl + (u), the eigenvalue 
of the transfer matrix acting on the reference state \<p) 

t{u)\4>) = S^X(u)A(u) + XX(u)4«(ti)} 10), (37) 



where we have introduced the transformations 

A» = A» - Sij^^-Aiu) = wr(«)A<(«)Ai(«), (38) 

with i — 1,2,3 and 

A(u) = uj a (u)A(u)A(u). (39) 



In the above expression, 



uj a {u) — 1, for m=l,...,5, 



(tt+2) («+£+) f j 

(u+l/2)(u+g++3) ' iUI 1 

(u+2)(uH+) f n _ / 

(u+l/2)(u+£ + +l) ' iU1 L 

(u+1/2) (-«+£+ -2)' 101 1 

(u+2)(u+U) f nr / 

(u+l/2)(«+£++2)> iUI 1 

u(-u+£--l) r _ 
(u+l/2)(u+£_)> iUI " L 

(ti+i/2)(-«-K_)' iU1 m 

u(-u+£_-l) f _ 
(u+l/2)(u+e_)' iU1 



= 1, 

= 2, 
= 3,4, 
= 5, 

1, 
3, 

2,4, 
5, 
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u 2 (u) 



(«+l/2)(u+C_)' 

tt(M+g-+l) 

(u+l/2)(-u+£_): 
(u+l/2)(-u+£_): 
I, u+1/2' 

M(-tt+g--l) 

(u+l/2)(u+£_)' 



^3 ( M ) = < 



u+1/2' 

u(ix+g_+l) 
(tt+l/2)(-u+£_)' 



TOT 
1UI 


III 


- 1 
— L i 


for 


771 


= 3 


for 


771 


= 4 


for 


771 


= 25 


for 


m 


= 1, 


for 


771 


= 3, 


for 


771 


= 4, 


for 


771 


= 2,5, 



I (u+l/2)(u+£-)' 

with ^(u) = #2 + (u), a#(u) = A"3+(u), c4(w) = fT4+(w) for / = 1, . . . , 5. In 
addition, 

A(u)A(u)|0) = w 1 (u) 2 \0), Ai(«)Ai(«)|0) = w 2 («) 2 |0). 

We note that the operators Bi(u), i = 1,2,3, constitute a three-component 
vector with both positive and negative pseudo-momenta still playing the role of 
the creation fields acting on the pseudovacuum state. In order to make further 
progress we return to the RE (10) and derive the commutation relations 



a ( \S / \ («2 - "1 + 1)(«1 + U 2 ) ~ , ~ ( 

A(«i)S (ii2) = — ■ 7 - T7 ———B a {u 2 )A{u 1 ) 



(u 2 - U\){ui + u 2 + 1) 
? ■ ——B a {ui)A{u 2 ) 

(U 2 - Ml) («1 + M 2 + 1) 



1 



Ml + U 2 + 1 

D bd (ui)B c (u 2 ) 



J2 B b(ui)D ba (u 2 ) + 5 ab 
b=\ 2u i + 1 

(ui — u 2 + l)(ui + u 2 + 2) 



S 6 («l)^(«2) 



(mi - m 2 )(mi + w 2 + 1) 
x { r «( Ul + « 2 + 1) Jr«( Ul - M2 )f d 4(«2)^(Mi)} 

4(tti + l)u 2 



(40) 



(41) 



(2«i + l)(2w 2 + + m 2 + 1) 



,(i) 



(2 Wl + l)SS ff ( Ul )i(« 2 ), 



between the diagonal fields and the creation fields. The summation convention 
is implied for repeated indices. The matrix which satisfies the Yang-Baxter 
equation, takes the form 



1, a = 1,2, 3, rf h 



1 



u+ 1 



, b = 1,2,3 
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f£ = — ,a^6 = l,2,3. (42) 

We notice that the first term in the rhs of each of the commutation relations 
(40)-(41) contribute to the eigenvalues of the transfer matrix, which should 
be analytic functions of the spectral parameter u. Consequently, the residues 
at singular points must vanish. This yields the Bethe ansatz equations, which 
in turn assure the cancellation of the unwanted terms in the eigenvalues of the 
transfer matrix after the whole nesting procedure. For convenience, we make 
a shift in the spectral parameters, u — v — 1/2, ui = Vi — 1/2, such that the 
eigenvalue of the transfer matrix (25) can be obtained as 



T(v)\$)=A(v,{Vi})\$) 

^- vn -v ^-;:^:^ m (43, 
- i/w,, - 1/2, n {v -;^t:; 1] ^. m>i»>. 

provided that 



W\{v t - l/2)WX(v t - 1/2) (2^ - 1) 
W?(vi - l/2)Wr(vi - l/2)(2v t + 1) 

^- (^-^ + l)(^ + ^ + l) A(1) 

^{Vi-Vi-lXVi + Vi-l) 

Ifr 

Here a = 1, 2, 3, and we appropriately choose 



WX(u) 



W+(u) = 



W~(u) 



1, for m — 1, ... ,5, 

(M+2)( " +g+) for I = 1, 

for I = 2, 

, for / = 3, 4, 



= < 



(u+l/2)(u+£++3) ' 

(u+2)(u+£+) 
(u+l/2)(u+£++l)' 

(u+2)(-u+£+) 
(u+l/2)(-u+e+-2) 

(u+2)(u+U) 
(u+l/2)(« 

1, for Z = 2,...,5, 

(«+i/5(u+$-)' for m = 1 > 



(45) 



(46) 



+S++2)' IOr f _ °' 

, for / = 1, 



u+l/2' 



for m — 2, ... ,5. 



(47) 
(48) 
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AW(i>, {vi}) is the eigenvalue of the second level transfer matrix t^> related 
to an SU (3) invariant open chain, i.e. 



r (1) =Tr K i +\v)T w {v)K i y{v)T {1 \v), (49) 



where 



TM(v, {v, i })=r$(v + Vi y h \l . . . rg\v + v N )ZZ~\ (50) 
fM(v, { Vl }) =r$}(v - v N )r N % • ..r$}(v - v 1 )^\ (51) 

Here we have used the standard notation where P is the 

standard permutation operator, which can be represented by a 3 2 x 3 2 matrix, 
i.e., p Q /3, 7 5 = S a s5/3 T It can be seen that the coefficients F l ,\{'"' ln act as the 
multi-particle vectors for the inhomogeneous transfer matrix (49). We remark 
that the coefficients W given in (45)- (48) are chosen in order to match the 
choice of the transfer matrix (49) with the nested K±^ -matrices, 



(Kl { l\v) 











K2 { l\v) 





V o 





K3<£ 



\ 



(52) 



Now corresponding to the first solution (20), we have 

Kl { l\v) =K2 ( l\v) = K3 ( l\v) = 1, (53) 
Kl { l\v)=K2 [ l\v) = K3 { +\v) = 1. (54) 

And to the second solution (21): 

K 1« W = K±M = 1, K2fl\v) = + 5- - 1 /. 2 . (55) 

v + £- - 1/2 

*!?>(,) = jra?>(,) = 1, K3?(v) = !;jl + !f /2 . (56) 
And to the third solution (22): 



/n«w=*2«w = ^ + *l + ^ 2 , K&iv) = 1, (57) 

A'l<>) = K2f{v) = 7_"/; + " 3 3 / 2 2 . = 1- (58) 
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And to the fourth solution (23): 



KlU(v) = 1, K2 { l\v) = K^\v) = V + i - + l {* (59) 



K\%\v) = 1, K2 { l\v) = K3$\v) = 7~/ + ,~ 3 / { 2 - (60) 
And to the last solution (24): 



-v + Z- + 1/2' 
-v-Z + - 3/2 
v - + 3/2 • 



*2«(t,) = l, Kl*\v)=K3*\v)= ^+£-W , (61) 

^(,) = 1, mg>(t,) = K#\v) = ";\ + _f /2 - (62) 

We can show that the reflection matrices (52) with the entries (53)-(62) do 
satisfy the RE 



r\ L 2 >(u - v) (u) r£(u + v) K ( _ L > (v) 



= K (1) ^ (v) r { S(u + v) K { } ] (u) r$(u - v), (63) 

1 tl 2 * 2 

r { S txt \v-u)K^ (u)r$(-u-v)KP (v) 

2 t 2 x ti 

= K| 1} (^rS^-u-u)/^ (u)rS )tlta (u-«). (64) 



3.2 Second-level Bethe ansatz 



In order to proceed in the nested algebraic Bethe ansatz, we have to repeat 
the whole procedure presented above for the internal block of the monodromy 
matrix. Similarly, we rewrite the second level transfer matrix (49) in the 
form 



T W{v) = Tr [kP(v)U { _\v)\ , (65) 

where U^\v) is defined by 

fW(v) = T^(v)K^{v)f^{v). (66) 

Now we label the elements of the monodromy matrices by 
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T«(A) 




BP 

u \\ 
r,(i) 

B? 

r,(i) 
-^11 
f,(i) 

-^21 

D (1) 
-^11 

r,(i) 

-^21 




(67) 



(68) 



(69) 



From the structure of the r (^-matrix (42), it is found that the operators B^(v) 
and Bjp(v) (a = 1,2) act as creation fields acting on the reference state 
1 



|0>* 



The operators Cf~\v 



, 2) behave as annihilation 



fields. In order to make further progress we return to the RE (63) and derive 
commutation relations, 



A {1 \ Vl )B^\v 2 ) = 
(vi ~v 2 )(v 1 + v 2 + 1) 



(v 1 -v 2 )(v 1 +v 2 + l) a KV yV 



B^( Vl )A^(v 2 ) 



(Ul - W 2 + l)(«i+U2 + 2) 



^1+^2 + 1 

x {r (2) (^i + v 2 + l)t h r^\ Vl - v^BW^D^v,)} 



2(«i + 1) 



.(2) 



(2ui + l)(ui-U2) 

4(ui + l)u 2 



(2^1 + i^^OAl 1 ^) 



(2ui + l)(2u2 + l)(ui + u 2 + l) 



+ 1 \2v 1 + l)%m\v 1 )AM(v 2 ), 



(70) 



(71) 



between the diagonal and the creation fields. Where again the summation 
convention is implied for the repeated indices. The matrix r^ 2 \v) is nothing 
but the SU(2) invariant i?-matrix, i.e. 



ba 



1, a 

V 

v + : 



1,2, 



ab 
ab 



V + l 



,a^b=l,2, 



a ^ b = 1,2. 



(72) 
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If we define the second level Bethe ansatz as 



(73) 



and make the rescalings Hi — > Hi — 1/2 and v — v — 1/2 we obtain from the 
commutation relations (70) and (71) the eigenvalue A^(v , {fi}{/ij}), i-e. 



= |<i(^ - 1/2)^ - 1/2) fi ^ - * - + * - 1} 



(v - /ii)(v + in) 



(74) 



wS(i?-l/2)WW(i;-l/2)n 



TV 



(6-^-l/2)(6 + ^-l/2) 



M 

*n 

Z=l 



- m + i)(v + in + 1) 



\ (v-Vi + l/2)(v + Vi + l/2) 

A (2) (^W,{M)ll^ (1) ), 



provided that 



W/j 1 ]^ - l/2jwgj(^ - l/2)(2/i, - 1) 
wjWi - l/^wlWi - V2)(2^ + 1) 

_-rr (/i/-^-l/2)(/i/ + ^-l/2) 

t -=i l/2)(/i/ + ^ + l/2) 



a; 



n 



(^ -/ij + ^(nt + Hi + 1) 



J (/!/ — At* — i)(m + fii - i) 



A^^W, {#}) \v= 



(75) 



Here a = 1,2, and 



_„ +f _ + l/ 2 > IOT "7, — <3, 

1, for m = 1,2,4, 

„+£_-l/2 ' IOr m _ °> 



(76) 



for Z 



f+3/2 
u+1/2' 

(i;+3/2)(tt+$ + -l/2) f , 

( v +l/2)( v +g + +l/2)' iUI f 

(^+3/2)(^+g + +l/2) f , 

(„+l/2) -3/2)' iU1 

(«+3/2)(-«+g++l/2) f , 



(j;+l/2) (-<;+£+ -3/2)' 
(«+3/2)(-«+g+-l/2) f j 
(v+l/2)(v+(++3/2) ' 101 4 



1, 

2, 
3, 
4, 
5, 



(77) 
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W$(v) = W$(v) = { 



(u+l/2)(t)+£_-l/2)' 

i) 

(tH-l/2)(-tH-£_+l/2) 

^+g_+3/2) 
(u+l/2)(-tH-£_+l/2) 

i> 

I. (t;+l/2)(t;+f_-l/2)' 

1, for I 

, for / = 3, 
4, 

for 1 = 5. 



for m 
, for m 
, for m 



for m 
■ 1, 



i_ 

-i>+£+-3/2' 
1 



I *>+£++3/2 ; 



2, 

3, 
4, 
5, 



(78) 



(79) 



Now A( 2 )(t>, {vi}, {/J-i}) is the eigenvalue of the third level transfer matrix 
related to an SU (2) invariant open chain, i.e. 



(2), 



(80) 



where 



T (2) (^, M, W)=r<8(v + ^)Z gi . . . r&\v + MET* ( 81 ) 
f {l \v, {vi}, W)=rjZ\v - HmYC^m ■ ■ ■ r^(v - (82) 



3.3 Third-level Bethe ansatz 



It has been shown so far that the eigenvalue problem of the transfer matrix 
is reduced to the diagonalization of the isotropic Heisenberg model with K- 
matrices 



v + Z--\ 


v + Z- o 

-v + Z 
1 
1 





v + S- 



, for 


m = 


2, 


for 


m = 


3,5, 


for 


m = 


1,4, 



(83) 
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((v + Z + 

-v + Z + , 
v + £+ + 1 

-v + £ + 

1 s 
1 



for 1 = 2, 
, for Z = 3, 5, 
for I = 1,4. 



(84) 



Following the derivation in [17], we immediately obtain the eigenvalues of the 
nested transfer matrix (80), given by 



-(2) 



Ii—Im 



a (2) M^MM,W)^ (1) 



Ii—Im 



(^-^)(^ + ^ + i) 



+^i 2 2 ) ^)^i 2 2 ) (^n 



Q 



1=1 \ 
provided that 



i=\ (v- m+ l)(v + fn + 1) 

(v - wi + l)(v + wi + 2) ) Fil) h-i M 
(v - wi)(v + Wi + 1) J 



(85) 



M 

n 



W% {wi)Wri (wi) (wi + 1) f=i ( w i - H + l)(w t + fij + 1) 



yr (Wj ~W m + l)(Wj + W m + 2) 
mil ( Wl -W m -l)(wi + W m ) 



(86) 



Here 



v + £_ — 1, for m = 2, 
■u + for m = 3, 5, 

1, for m = 1,4, 

for 1 = 3,5, 
for Z = 1,4, 



(87) 



2S+1 

v±l 
6+1/2' 



(88) 



5+1/2 > 101 — Z > 

-g(g-g-+l) 



i5+l/2 



g^.for m = 3,5, 
for m = 1,4, 



(89) 
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- 2, for 


1 = 2, 








— 1, for 


I = 3,5, 


(90) 




1, 


for 


I = 1,4. 





The eigenvalues (74) and (85) as well as the constraints (75) (86) on the 
rapidities jii and wi have paved the way for the complete diagonalization of the 
transfer matrix (25). Making a further shift on the rapidities, Wi — + Wi — 1/2, 
w m — > u> m — 1/2 and {; = v + 1/2, the eigenvalues of the transfer matrix (25) 
are given by 



A(v, { Vi }, M, { Wj }) = Kl.(v - \)Kl + {v - l -) 
{Wt(v -\)WX(v -\){v + \f L 

N 

><n 



1=1 



(V-Vi)(v + Vi) 



+ W+(v - \)Wr(v - \)W$(v)w!$(v)(v - \) 2L 

" ( V - Vi + l)(y + Vi + j) M ( v _ U[ _l^ v + Ui _ 1) 
X L\ ( V - V i)( V + V i) t=l (V-Ut + \){V + Ui + d 
+W+(v - l)Wf(v - \)W${v)wV(v)W$(v + \)W™{v + \){v - \f L 

x yr (V^JM+ j)(v + IM + f) (v - Wj)(v + Wj) 

1=1 + \){V +m + \) fJi (V - Wj + l)(v + Wi + 1) 

+W+(v - \)W 3 -(v - \)W$(v)wV(v)W®(v + \)W${v + \){v - \f L 

x ^( r w 3 + 2)(v + w J + 2)\ 
i=i (v- Wj + l)(v + wi + l) [ V ; 



The three rapidities {vi, /ij, Wk} of flavor waves satisfy the Bethe ansatz equa- 
tions 



1\2L 
2> 



[Vi - 



l\2L 



N 

n 



(Vi- Vl + 1)^ + ^ + 1) g( 



5)(«i + - \] 



,tl («i -Ui - + - 1) , =1 

# 1 ,w* i ,c-)n ( " i ""' +l)fe+ " ,+l) 



+ 5)(ui + A*i + 



1 )' 

2^ 



l= l 



\)(Hj + Vi - |) 



(92) 
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il (flj - fjq + VjjfjLj + flj + I) ® (flj -Wi- +W t -\) 
11 (,.. _ ,,. _ IV,,. _L_ ,,. _ 1\ 11 



hi + |)(^fe + hi + 1) 



(93) 



i=i k 



n 



Hi - \){w k + Hi~\) 
w t + l)(w k + w t + l) 



\ (w k -wt - l)(w k + wi-iy 



(94) 



i^k 



for % = 1, . . . , TV, j = 1, . . . , M and k 
introduced the notation 



, Q, respectively. Here, we have 



C(«i,£fc) 



v(Hj,^ 



( ^i±|±^i, for Z = 1, 



t, for I = 3, 



m 
m 



1, 

3, 



1, for Z = 2,4, m = 2,4, 



f, for / = 5, m = 5, 

for Z = 1, 2, 3, m = 
f, for / = 4, m = 4, 



1,2,3 



for I = 5, m = 5, 



, for / 



for J = 3,5, 
1, for Z = 1,4, 



m = 2, 
m = 3, 5, 
m 



1.4. 



(95) 



(96) 



(97) 



These boundary factors coupled to the three degrees of freedom will result in 
a rich physical scenario. From (27) and (91), we finally obtain the eigenvalues 
of the Hamiltonian (1) as 



E = < 



+ J2l 1 (^^ L T + J±), for l = m = 1,2,5, 

r ,Jw 



(98) 



N 



i(^^I 4 +J±), for Z = m = 3,4. 
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4 Boundary impurity effects 



Having diagonalised the Hamiltonian (1) by means of the algebraic Bethe 
ansatz, the next step is to derive the thermodynamic Bethe ansatz equations. 



4.1 Derivation of TBA 



For later convenience in the analysis of the Bethe ansatz equations, we make 
the change of variables: Vi — > —ivi, fii — > —i/ii, Wk — > —vwk and some rescalings 
in the boundary parameters £±. The Bethe ansatz equations are now 



N . M i 

Vi - rvi - 1 yj Vi - ru t + 5 



i \2L 



t( r \ TT TT Yi_!_zi^l± TT _ l: I 



r=± ;=i 
M 



N 



■kw. 1111 .11 - + f n w ~ "" + 5 



r=±i 



± f=I ^' ~ r ^ + 1 i=i ^' ~~ rw i ~ h 7=1 H ~ rVi ~ h 

Q . M 



r =± r=i ~ ru,/ + 1 t=i Wk ~ r ^ i ~h 

l^k 



(99) 



(100) 



(101) 



where 



^±i|±, for Z = 1,3,5, m = 1,3,5, 

Vi—ip± ' 

1, for Z = 2,4, m = 2,4, 



1, 

Itj-i0± ' 



for / = 1,2,3, m = 1,2,3, 
for I — 4, m = 4, 



n(w k ,p±) 



L Srl§' for / = 5 ' m = 5 ' 

5±g, for 1 = 2,5, m = 2,5, 
^Ti^+li , for / = 3, m = 3, 



1. 



for / = 1,4, m = 1,4. 



(102) 



(103) 



(104) 



The shifts in the parameters £± are given by 
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l3± = i±-\, for/ = m = 1,4,5, 

/3± = £t-§, for/ = m = 2, (105) 
/3± = — £± — \, for / = m = 3. 

Correspondingly, the energy spectrum is given by 

^t(-7?V4 (106) 

Here we have dropped some constants appearing in Eq. (98), which will be 
used in deriving one point correlation functions later. 

From the above Bethe ansatz equations (99)-(101), it is found that in the cases 
I — m — 1,3,5, the solutions v\ = ±i/?_ and v T = ±i/3 + form two boundary 
bound sates in the charge rapidity when (3± are negative. Nevertheless, in the 
case / = m = 5, besides the charge boundary bound states, the boundary 
bound states exist also in the spin rapidites, i.e. 



(±i(/3- + §), (±i/3_, 
(±i(/3+ + |), [±i/3+. 

No boundary bound state exists in the remaining cases. We observe that when 
J± > — - (1 — cosk), the reference state becomes the true ground state, i.e., 
the ground state is given by a product of the singlet rung states. The minimal 
gap can be easily calculated and is given by 

A = J ± - ^(1 -cos A;), (107) 



where k = 7r /[l + 4j;(^ + ^ _ )] - It is obvious that gap remains almost unchanged 
in the thermodynamic limit and is almost the same as A = Jj_ — in the 

periodic case because L >> In the regime — | < (3± < ~\\jl — ^7^, the 

boundary bound states are stable. Otherwise, in the remaining regime, they 
become excited states. In the limit J± — > oo, all the boundary bound states 
are excitations. We shall see that the boundary bound states radically affect 
the edge ground state properties. For J~ = — i(-^= — In 3) < J± < the 
ground state consists of three branches of Luttinger liquids associated with 
the rapidities v, fi and w. Here J~ is the critical transition point from the 
577(3) phase into the 577(4) phase in the absence of a magnetic field. The 
triplet states can exist in the ground state. This corresponds to a continuum 
of massless excitations. 
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The thermodynamics of the boundary fields can be derived from the Bethe 
ansatz equations (99)-(101). We now focus on the analysis of the Bethe ansatz 
equations. As usual, we define the functions 



x + i- 

e n (x) = r|, n (x) = ilne n (x), 



In Id 

a »M = S?Tf s SS 9 " (I) ' (108) 

in terms of which the Bethe ansatz equations (99)-(101) become 



M N 

av t ,p±)(e 1 (v l )) 2L n iw^-^) = n iw^-™*), 

r=±Z=l r=±j=i 
AT M Q 

v(»j,p±) im^ifij -wi) = n ri e 2(^ - r ^)ii e -i(^ - r ^)( n o) 

r=±l=l r =±i=i 1=1 

1+3 

M Q 

Sl{w k ,P±) II IT ei(w fc - rm) = Yl Y[e2(w k - rwi). (Ill) 

r=± 1=1 r=± i = i 

l^k 

In order to study the thermodynamics of the model with boundary impurities 
we begin by adopting the string hypothesis [26]. If we define V-j = —Vj, 
fi-i = —fj,i and W-k = w k , the Bethe ansatz equations (99)-(101) admit the 
string solutions 



v n aiJ =v n ai +i-(n+l-2j), 

in thermodynamic limit, where j = 1, • • • , n, a a = 1, • • • , and v™ v u™ 2 
and w™ 3 are the positions of the center of the strings. The number of n-strings 
Njjf> satisfy the relation = J2 n nN^ a \ By taking the thermodynamic limit, 
the Bethe ansatz equations become 



ti )h = a n + - E Anm * P $ + E ^ * pS , ("2) 

m m 

P ( : )h = - E Anm * pS? + E «™ * (P« + P®), (H3) 
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PI 



(3)h _ J_ (3) _ A 



2L 



m 



* p 



(2) 



(114) 



where the symbol * denotes the usual convolution. Here p$(v ), a = 1, 2, 3 are 
the densities of roots of the three flavors, p^ h (v ), a = 1, 2, 3 are the densities 
of holes of the three flavors and p^ n , i = 1,2,3 are the contributions from 
boundary fields associated with different rapidities. These boundary phase 
factors are given by 



o (l) 

Pbn 



= < 



E±Er=i a r i +2 / 3±+i-2«(^) + a„ 2 (A), for 1 = 1,3,5, m = 1,3,5, 



(115) 



o {2) 

Pbn 



(3) 
Pbn 



In addition 



(116) 



a n2 (A), for I = 2,4, m = 2,4, 

a n2 (A) - a„i(A), for / = 1, 2, 3, m = 1, 2, 3, 

E± £p=i an+2/9±+i-2j(A) + a n2 (A) - a„i(A), for / = 4, m = 4, 

- E± E™=1 «ra+2/3 ± + 2 - 2 «(A) 

+a n2 (A) - a„i(A), for 1 = 5, m — 5, 

E±Er=i a n+2/3 ± +i-2«(A) + a n2 (A), for Z = 2, 5, m = 2,5, 

- E± Er=i On+2 i 9i+3-2j(A) + a n2 (A), for / = 3, m = 3, 
a n2 (A), for 1 = 1,4, m = l,4. 



117) 



Aim(A) = 5(A)5„ m + (1 — 5 nm )a|„_ m |(A) + a„ +m (A) 

Min(n,m)-1 

+2 ^2 a|„_ m | +2i (A), 
i=i 

Min(n,m) 

(A)= E (A). 

We emphasize that the boundary potentials enter in the expression for the 
ground state energy implicitly via p^ (v) in the above equations, with contri- 
butions to the densities of the roots at the order of 1/L. In order to find the 
equilibrium state of the system at fixed temperature T and external magnetic 
field H (> 0), we minimize the free energy F = E — TS — HS Z with respect 
to the densities to obtain the TBA in the form 



Hl + Vn 2) ) 

\ln(l+J?i 3) ) J 



T 



Em Anrn ~ Em a ™ m ^ 
~~ Em a nm Em ^nm ~ Em a nn 
\ — Em a nm Em ^nm J 



Vm 

Mi + -w) 

^Ml + 4y)J 



.(118) 
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The driving matrix G n depends on the choice of the reference state. Explicitly, 
for Jj_ < 0, G = column(— '^-2ira n + nH,nH, —n(J± + H)), giving the free 
energy 



F(T,H) 
2L 



/ ^ MA) 
^a n (A) ln(l + e— 

-oo n=1 

rp 3 °° OO (a) 

^E / EpS ln(l + e-^)dA. (119) 



a=l_„ n=l 



On the other hand, for Jj_ > 0, G = colum(— ^-2na n + n(J± — H),nH,nH) 
and the free energy is given by 



£eia = _ r /f a . (A)M1+e -^ )<tt 

^ J n=l 

-oo " 1 

-2lE / EpSm + e-^)dA. (120) 

a=l_oo n=l 

Here ^(A) = p® h (\)/p®(\) = exp(e^(A)/T), / = 1,2,3, with the dressed 
energy e$ playing the role of an excitation energy measured from the Fermi 
level. 



Using the relations, 



(oo + a 2 ) * InijW = ai * [ln(l + £1) + In(l + V^i) 

-Ml + ^-lna + ^pry), (121) 

Tjn Tjn 

another form of the TBA is given by 



e (<0 00 £ < a ) 

ei a) = 9i ] + Ta 2 * ln(l + e"^) + T(a + a 2 ) E a ™ * ln ( X + e"^) 

m=l 

OO / e (o-l) e (o+l) \ 

— T E a ™ * ( + e-^) + ln(l + e~^) , (122) 



m=l 



>) 
e n-l 



ei a) =3i a) + Ta 1 *ln(l + e^) 

(a) OO (a) 

+Ta 2 * ln(l + e - V) + T(a + a 2 ) E a ™-« * ln ( 1 + e ~^) 

m>n 
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OO / e (o-l) £ (a+l) \ 

—T £ a m _ n+1 * ln(l + e"-^) + ln(l + e"^) , (123) 

m>n \ / 

for n > 2. In the above a = 1,2,3 and e^(A) = e^(A) = is assumed. The 
driving terms are given explicitly by 





= _!±2na 1 + H, 


= h, 








9? 


= H, 


# = H, 


for 


J± 


<0, 


9? 


= -(J ± + H), 


g® = -(J± + H), 








»!■' = 


- J -±2n ai + J ± - 


H, gW = J ± -H, 








(2) 

9i = 


H, 


# = H, 


for 


J± 


> o. 


(3) 
#1 = 


H, 


# = h. 









(124) 



(125) 



4-2 Boundary bound states and impurity effects 



In the low temperature limit, the states with positive dressed energy are empty. 
The zeros of the dressed energies define the fermi energies. We decompose 
into its positive and negative parts, = e^ + +e^~ . An analysis of equations 
(122) and (123) in the limit T — > reveals that for the ground state, the roots 
are all real corresponding to n — 1. All dressed energies with n > 2 
correspond to excitations. Thus the TBA for the ground state is, for J± < 0, 



g?-a 2 *e^- 



+ a\ * e 



(2)- 



e ( 2 ) = H -a 2 * e (2) ~ + a x * 



+ e 



(3)- 



e (3) = _#_ J ± _ a2 * e (3)- +ai * e (2)- ; 

and for J± > 0, 



(126) 



9 ? - a 2 * e«- 



e 

e (3) =J fJ_ a2>Ke (3)- 



ai * e 



(2)- 



+ a 1 * 
+ a x * e 



(i)- 



r(3)" 



(127) 



In this case, the free energy is given by 



F(0,H) 



2L 



-H + f? 0O a 1 (\)e?-(\)d\ + ±f h , for J ± < 0, 



roo«i(A)4 lK (A)rfA + ^/ b , 



128) 



for J ± > 0, 



26 



where 



/b = £ / pS(A) e S aK (A) + 0(/? ± + p e )Eh 

a=1 -Qa 



Qa 



(129) 



and 9(x) denotes a step-like function. Define (3 C 




— r- then in the 



interval — \ < (3 < (3 C) 9(f3) = 1, else 9(f3) = 0. In the above E hs denotes the 
boundary bound state energy, given by 



It is worth noticing that if f3± < f3 c we should take the boundary bound states 
into account in the boundary contributions p^} for the cases I — m — 1,3,5. 
The TBA (126) and (127) provide a clear physical picture of the ground- 
state and in turn the thermodynamic properties, such as the free energy, the 
magnetization, the susceptibility, etc. The boundary impurities coupled to 
the three rapidities affect the low temperature physics at the edges in various 
different ways, which we now explore. 

From the TBA (127), we notice that if J± > = the triplet excitations 

are massive with energy gap A = J± — The rescaling 7 = 4 was fixed [14] 
for strong coupling compounds, e.g. (^(CsH^^^CL, [3], (CsH^N^CuB^ 
[4], etc. Here J+ is the critical point at which the quantum phase transition 
from the three branches of Luttinger liquid to the dimerized U (1) phase occurs. 
If J_|_ > J+, we can show that in the presence of a strong magnetic field two 
of the triplet states (|3) and |4) in (5)) in the bulk part will never be involved 
in the ground state. However, at the boundaries this is not always true due 
to the presence of the boundary impurities. In a strong magnetic field the 
ground-state may be considered as a condensate of SU (2) hard-core bosons. 
The gap is reduced by the magnetic field H, i.e. A = Jj_ — — H. Thus 
the first critical field occurs at the point H ci where the gap is closed, i.e. 
gfiB^d = J± — The quantum phase transition from a gapped phase to 
gapless Luttinger phase occurs. By continuing to increase the magnetic field 
H over H cl , the triplet state |2) becomes involved in the ground state with 
a finite susceptibility, also affected by the boundary impurities in the low 
concentration regime. If the magnetic field is greater than the rung coupling, 
i.e. h > J±, the triplet component |2) becomes the lowest level. Therefore, it 
is reasonable to choose the basis order as (|2), |1), |3), |4)) T . Subsequently the 
driving terms are given by = — 2irJ\\a 1 — J±+H, = J± and = H. A 
second critical field H c2 (H c2 > H c i) can be determined by the magnetization 




for Jl > 



for Jl < 0. 



(130) 
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arriving at its saturation value S z = 1. Then the reference state becomes the 
true physical state and the critical field H c2 is given by 

H c2 = J ± + ^. (131) 

7 



In this case, all the boundary impurities are gapfull with the ferromagnetic 
gap A = j2 B g(H - H c2 ). 

Let us now discuss the boundary impurity effects in the vicinity of the critical 
point H c \. After a lengthy calculation, similar to that employed in [14] for the 
periodic case, we find the free energy in the presence of a strong magnetic field 
H, 

F(0,H) „ 4Q(J+ - J cff ) / _ 2Q\ 1 

2L 7T V 7T J 2L 



(132) 



where Q is the fermi point given by Q ~ y 4 j+_ c 5 ( jJ l ff g ^ an d /b is the surface 
free energy from the boundary impurities in the vicinity of H c \. Explicitly, 
for \j3±\ > |, or say < U± < or U± < 0, it is given by 

_ 2Q(j+-j eff ) / 1 + i + i n fo r j = m = i 3 5 

/b~{ * ^ + ^ (133) 

■ _ 2 Q( ^-j off ) > for / = m = 2?4 . 



For -i < /3 ± << /3 C , or say C/± >> C/ bs = 2JJj(l - Jl - -^-), we have 



7J ± ' 



/b 



_ V 2Q(J+-J eff ) / 3 , 1 , 1 1_\ 

^± w \2 /3± /3±+l /3±-lJ 

+ E ± (-ir^ + i), 

2Q(J+-J eff ) 



for / = m = 1, 3, 5, (134) 
for Z = m — 2, 4. 



While for |/3±| very small, or say ^ < U± < U^ s , it is given by 




+ £±/(#fc) for / = m = 1,3,5, 

(135) 



where 
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8J11 



7T 7 1 - 4pl 



wcta « 4Q/J± + ^I arcta « 

, P± * P± 



4 ^ + -^ ff ) 1 / Q_ _ \ 



(136) 



In the above J c s = J± — H and the parameters f3± are related to the boundary 
impurity coupling U± by 



1 



2U± 



— , for/ = m = 1,4, 5, 



7 

2[/± 



2U ± 
l -+U ± '' 



for / = m = 2, 



for / = m = 3. 



(137) 



The magnetic susceptibility follows from \ ~ ~l^^?iP~' Here, to illustrate 
the boundary effects, we will focus on the discussion of the strong coupling 
compounds J± » J\\ with the boundary impurities in the case / = m = 1. 
Other regimes can be handled in a similar way. It is very clear that the stable 
boundary bound states are exhibited only in the strong ferromagnetic bound- 
ary coupling U± > Ubs- In Eq.(137), we emphasize that the mathematical 
singular points do not exist, or alternatively f3± — does not mean that the 
rhs of Eq. (137) has singular points. For instance, if U± = 2J\\/j, the boundary 
parameters £± = 1/2. Thus the phase factors in the Bethe ansatz equations 
(99), (100) and (101) are equal to 1 for the case / = m — 1. In such a case, 
the model exihibits special symmetry (the quantum algebra SU q (4) invariant 
Bethe ansatz equations) which leads to a different expression for the bound- 
ary free energy than the above ones. For antiferromagnetic boundary coupling 
U± < 0, the susceptibility is given by 




X* T— = = 1 + 77 E(o + ~) . ( 138 ) 
ny/AJ+(H - H cl ) 



while for the strong ferromagnetic coupling U± » Ub s , with U± > 0, 



Notice, in both cases, that the susceptibility diverges with the square root of 
the field in the bulk and in the boundaries. In addition, the susceptibility at 
the boundaries is enhanced or decreased by different impurity dopings. This 
behaviour is illustrated in figure 2. From the Bethe ansatz equations, we can 
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Fig. 2. The susceptibility versus magnetic H for different impurity coupling U± = U 
(U = corresponds to free boundary conditions). Here we consider the strong 
coupling compound Cu2(C5Hi2N2)2Ci4 [3] with J± = 13. 2K, Ju = 2.5K and 7 = 4 
with an impurity concentration 2 percent in a ladder with length L = 50. 

also calculate the magnetization in the vicinity of H c i, 



2L 



7T 7T IL ± 



(140) 



For antiferromagnetic boundary coupling U± < this expression reduces to 

(141) 



S .. E !9 (1 4)(LL 



± 



7T 



± 



while for strong ferromagnetic boundary coupling U± » U^ s with U± > 

h/ — d-rrHo I" i — r-- J 



(142) 



A plot of the magnetization S z against the magnetic field for different bound- 
ary impurities U± is given in figure 3. By analyzing both figures we can observe 
the competition between the boundary impurities and the magnetic field in 
the thermodynamic properties. In particular, we find an enhancement of the 
susceptibility in the weak anti- and ferromagnetic regimes (we consider the 
sizes 2L = 100, the impurity concentration 2 percent ). The susceptibility and 
the magnetization are lifted slightly in the weak antiferromagnetic bound- 
ary regime in the case of open boundary conditions, whereas they contribute 
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0.5 



bulk part 




10.5 10.7 10.9 11.1 11.3 11.5 11.7 11.9 

magnetic field H 

Fig. 3. The magnetization versus magnetic field H for different impurity coupling 
U± = U (U = corresponds to free boundary conditions). As in figure 2, we consider 
the strong coupling compound Cu2(C5Hi2N2)2Cl4 [3] with J± = 13. 2K, Jm = 2.5K 
and 7 = 4 with an impurity concentration 2 percent in a ladder with length L = 50. 

negatively to the bulk when U± becomes more and more negative. This is 
reasonable, since negative U± energetically favours the singlet state (recall the 
boundary terms in the Hamiltonian (1)), even if the magnetic field is very 
strong, such that the spin-1 component of the triplet is involved in the ground 
state. The point is that a very negative U± can overcome the spin-1 component 
of the triplet and dominate the edge rung state. In this circumstance, the edge 
state is a pure singlet state and the edge magnetization (141) is zero due to the 
fact that U± effectively decreases the edge magnetic field H to H c i such that 
the fermi boundary Q = 0. This results in negative susceptibility and magne- 
tization contributions to the bulk. In contrast to this case, the ferromagnetic 
impurities lift the susceptibility and the magnetization in the weak coupling 
regime U± < When U± becomes larger, the triplet edge state is energeti- 
cally favoured so that the boundary coupling can overcome the magnetic field 
to bring the three components of the triplet into the edge states. Therefore it 
causes a negative contribution to the bulk susceptibility and magnetization. 
This situation is different from the case of the bulk impurities, where the sus- 
ceptibility is increased by the impurity coupling due to the forward-scattering. 
This fact can be seen clearly from the one point correlation function of the 
ground state at the edges, for antiferromagnetic boundary coupling and weak 
ferromagnetic boundary coupling, i.e. U± < — , 



{Sa ' Ta) + dU^ fh + * (2J^-U ± r ( } 

and for ferromagnetic impurities in the strong coupling regime U± » Pbs, 
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Fig. 4. One point correlation function (143) versus antiferromagnetic boundary cou- 
pling U as a function of magnetic field. The curve is lifted by the magnetic held, how- 
ever it is decreased by the boundary impurities which favour the singlet state. Here 
we consider the strong coupling compound Cu2(C5Hi2N2)2Cl4 [3] with J± = 13.2K, 
J|| = 2.5K and 7 = 4 and U± = U. The case U = corresponds to the free boundary 
effect. 



(S a .f a ) 



1 

4 



2Q(H - H, 



cl 



71 



4J, 



4J||/ 7 



(2J|, 



3fA 



+ 



(2J||/ 7 -C/±; 
4J||/ 7 



(2J,|/ 7 + C4 



:i- 



J\\ ■ 



(144) 



In the above a — 1, L. The boundary one point correlation functions are given 
by 



(S a .f a ) = -^(N S )+ 1 -(N T ). 



(145) 



Here N s and iV T are the probabilities of the singlet and triplet state respec- 
tively. This is because the eigenvalue of the one point correlation function 
(S a .T a ) acting on the singlet (triplet) state is — | (■?). We have plotted the 
correlation function for antiferromagnetic boundary coupling in figure 4. 



We see that the magnetic field always lifts the spin-1 triplet component. How- 
ever, in the case of antiferromagnetic boundary impurities and open bound- 
aries the singlet state is favoured as long as U± becomes more negative, the 
triplet moves out of the edge state and the one point correlation function tends 
to — |. On the other hand, for ferromagnetic coupling impurities (see figure 5), 
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strong ferromagnetic coupling U 



Fig. 5. One point correlation function vs ferromegnetic boundary coupling U :(a) The 
function (143) is lifted by the magnetic field and weak magnetic impurity coupling 
U. (b) The function (144) tends to | as the boundary impurity coupling becomes 
larger. Here we again consider the strong coupling compound (^(CsH^^^CLi [3] 
with Jj_ = 13. 2K, J|| = 2.5K and 7 = 4 and U± = U. The case U = corresponds 
to free boundaries. 



the correlation function increases due to the ferromagnetic properties and the 
magnetic field in the weak coupling regime U± < However, if U± becomes 
larger, the three components of the triplet get involved in the edge state, such 
that the correlation function tends to \ for strong ferromagneitc impurity 
coupling. This result indicates that the edge state can be a pure singlet state 
in the strong antiferromagnetic boundary coupling regime whereas it turns 
out to be a pure triplet state in the strong ferromagnetic boundary coupling 
regime. This reveals the role of anti- and ferromagnetic impurities. 



On the other hand, the boundary impurities coupled to the spin degrees of 
freedom, namely, and p h will also affect the ground state properties non- 
trivially. From the free energy (128), these impurity densities will contribute 
to the low energy. Considering the case J± < 0, in the absence of the magnetic 
field, the triplet is completely degenerate while the fermi surface of the singlet 
is lifted as J± becomes more negative. Certainly, if J± < J~ = — ^(75 — In 3) 

the singlet rung state is not involved in the ground-state, namely e^ 3 -*(0) > 0, 
whereas two triplet fermi seas still have fermi boundaries at infinity. Under 
such a configuration, the dressed energy potentials are 



(A) 



27rJ|| cosh I A 
a/37 cosn 71 A ' 



(A) 



2ttJ|| sinhfA 
a/37 sinh ttX 



(146) 



33 



The free energy can be given by 

^,-> 1 )-,(i,) + ^, 



where 



h= j pff(A) e i 1) (A)dA+ j f$(\)e?\\)d\. (148) 

— oo — oo 

The first part in (147) is nothing but the standard £77(3) ground state energy 
of the bulk. The remaining part is the boundary surface energy for various 
boundary impurities. 



5 Conclusion and discussion 



In summary, we have discussed in detail the algebraic Bethe-ansatz solution of 
an integrable spin ladder system based on the SU (4) symmetry with bound- 
ary impurities. Five different classes of solutions of the graded RE leading 
to different boundary rung interactions in the Hamiltonian were obtained. 
The Bethe-ansatz equations, the eigenvalues of the transfer matrix and the 
energy spectrum were given explicitly. Furthermore, the three-level transfer 
matrices, characterizing the different flavour sectors separately, allowed us to 
embed different impurities into the system. From the Bethe ansatz solutions 
(112)-(114), we found that the boundary impurity effects characterized by 
C(vi,£±), v(fJ>j,£±) and fi(iUfc,£±) act indeed non trivially on the densities of 
roots for the three rapidities and thus change the ground state properties, 
the boundary bound states as well as the low-lying energy spectrum. In the 
thermodynamic limit, the spin gap remains almost unchanged. However, the 
boundary susceptitblity and magnetization reveal novel magnetic properties 
for strong and weak impurity couplings. In strong impurity coupling, the im- 
purities induced by the open boundary conditions can result in either a pure 
triplet or a singlet edge state due to the nature of the pure back-scattering at 
the edges and the magnetic impurities. Strictly speaking, the edge state can 
be a pure singlet state in a strong antiferromagnetic boundary coupling regime 
whereas a triplet state with an effective magnetic moment can exist in a strong 
ferromagnetic boundary coupling regime. Correspondingly, the one point cor- 
relation function for strong antiferromagnetic boundary impurities tends to 
the singlet eigenvalue — |, whereas for strong ferromagneitc impurity coupling 
it tends to the triplet eigenvalue \. This behaviour may be observed in ex- 
periments due to different boundary magnetic moments. Although the TBA 
solution of the SU (4) ladder model (1) predicts the quantum phase diagram in 
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good agreement with experimental results for the strong coupling compounds, 
the full finite temperature thermodynamic properties of the model remain to 
be calculated. 
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